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Radio observations of the nova remnant of GK Per around 1985 (epoch 1) have shown it to have a spectral 
index of ^ —0.4 {S oc v"') at frequencies < 1.4 GHz and a spectral index of —0.7 at frequencies > 1.4 GHz 
CJ, |5 |. This spectrum is either due to secondary electrons from p-p interactions |5 1 or due to the presence of 
at least two electron populations which dominate the emission at different frequencies. Our radio observations 
(21 of the remnant of GK Persei using the Giant Metrewave Radio Telescope (GMRT) and archival Very Large 
Array (VLA) data show that the nova remnant has a low frequency (< 1.3 GHz) spectral index a ^ —0.85 
(in 2003; epoch 3) and a high frequency spectral index of a ~ —0.7 (1997; epoch 2). Moreover the higher 
frequencies show an annual secular decrease of ^ 2.1% between epochs 1 and 2 whereas there is no change 
at 0.33 GHz. This suggests a significant evolution of the spectrum over 20 years. In this paper, we discuss a 
possible physical basis for this spectral evolution based on a comparison of the synchrotron loss and electron 
acceleration timescales. 

1. Introduction 

GK Persei is a classical nova system that had a nova outburst in 1901 . Classical novae are a type of cataclysmic 
variable objects which are characterized by a significant mass loss (~ 10^^ — IO^^Mq) during an outburst 
which are infrequent (once in about 10^ — 10^ years). GK Per consists of a white dwarf primary accreting 
matter from an evolved K2 subgiant secondary 1 3 1 and the system is embedded in an extended bipolar nebula 
Q> probably a leftover from its earUer avataar and which continues to influence the evolution of the nova 
remnant from the 1901 outburst. 

The nova remnant has been extensively studied in the optical f2l (and references therein). Radio observations 
at 1 .5 GHz and 4.9 GHz by |6J discovered the presence of a non-thermal, polarized radio remnant, with a 
morphology resembling the optical remnant. A detailed multifrequency study of the nova remnant and the 
environs of GK Per were made by fT] which included radio continuum study at frequencies between 0. 15 1 and 
4.86 GHz, in which they found that the radio spectrum had a low frequency turnover around 1 GHz, with the 
spectrum being flat at the lower frequencies. Biermann et al. |5 1 have reanalysed some of the Seaquist et al. Q 
data and find that the low frequency spectrum has ^ —0.4. They suggest that the origin of this spectrum lies 
in the p-p interactions. The radio emission and the optical images of the remnant indicate that the nova ejecta 
are interacting primarily with the ambient medium located to the southwest of the nova. The optical nebular 
remnant shows a deceleration in this direction. 

Anupama & Kantharia (2005) recently studied the nova remnant in the optical using the 2-m Himalayan Chan- 
dra Telescope (HCT) and in the radio continuum at 0.33, 0.61 and 1.28 GHz using the GMRT (2003; epoch 3). 
They also used VLA archival data at 1.49 and 4.86 GHz obtained in 1997 (epoch 2). Anupama & Kantharia 
I2I find an annual secular decrease of about 2.1% in the flux density at 1.49 and 4.86 GHz. They also find 
that there is a break in the spectrum near 1 GHz with the spectrum getting steeper at the lower frequencies, in 
conflict with the Seaquist et al j7| (epoch 1) result. Here, we try to explain the results of Anupama & Kantharia 
(2005) and comment on the contribution of GK Per-like novae to the total cosmic ray content. 
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2. The radio emission 

Fig. ^shows the GMRT images superposed on the optical [Nil] HCT image. The optical and radio shells are 
of similar extent. The low-brightness radio emission extends to the north-east where Anupama & Kantharia 
(2005) report enhanced density. The peak emission at 0.33 and 0.61 GHz appears to arise from a location that 
is different from that of the peak emission at 1 .4 GHz. This might be indicating the presence of at least two 
dominant electron populations in the nova remnant or different physical processes dominating at the different 
frequencies. The spectrum of the high frequency peak is slightly flatter compared to the low frequency peak 
probably indicating local acceleration. 




Figure 1. GMRT radio images (contours) of GK Per nova remnant shell at 0.33 GHz (left) and 0.61 GHz (middle) and 
VLA 1.4 GHz (right) superposed on the [Nil] optical image (grey scale) from HCT (From [2 J). 

The radio spectrum of the nebular remnant of GK Per which includes data from Seaquist et al. Q and Biermann 
et al. who revised the flux estimates given by Seaquist et al. at 0.33 and 0.61 GHz is presented in Fig.|2ja). 
For convenience we refer to the data from Seaquist et al. and Biermann et al. taken around 1985 as epoch 1, 
the archival VLA data from 1997 as epoch 2 and our data in 2003 as epoch 3 in this paper From the figure, it 
is clearly seen that the flux at all frequencies has decreased with time, except at 0.33 GHz, where the flux has 
remained unchanged. The decrease in flux density at 1.49 and 4.8 GHz can be explained by a secular decrease 
of 2.1% per year [2] which has left the spectral index unchanged at —0.7 whereas at 0.6 GHz, the decrease is 
about 1.7% per year. The object which had a spectral index of —0.7 between 1.4 and 4.8 GHz and a flatter 
index of —0.4 at frequencies lower than 1.4 GHz (Seaquist et al. 1989, Biermann et al 1996) seems to have 
evolved in the last 18 years giving a dramatic steepening of the low frequency spectrum to an index of -0.85 for 
frequencies lower than 1.3 GHz (21 These observations raise questions such as what could be the underlying 
physical reason for the two spectra which we discuss in the next section. 

3. Spectral evolution and electron acceleration 

In epoch 1, the spectrum had a spectral index —0.7 at frequencies above 1.4 GHz and an index of —0.4 at lower 
frequencies. The most convincing explanation we could find in the literature was due to Biermann et al. (5). 
Assuming that the electrons responsible for the radio emission were secondary products of p-p collisions, they 
derived a high frequency spectral index of —0.7 and a low frequency index of +1/3. It may be noted that the 
low frequency spectral index that Biermann et al. |5J predict does not match the observed value. 

As mentioned earlier, we find that the spectrum above 1.4 GHz has undergone a secular decrease between 
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epoch 1 and epoch 2. There are no measurements below 1.4 GHz during epoch 2 and no measurements above 
1.4 GHz during epoch 3. We find that the low frequency spectrum (below 1.3 GHz) has a considerably steeper 
spectral index (—0.85) during epoch 3. Assuming that the high frequency spectrum continues to follow the 
same secular decrease observed between epochs 1 and 2, we have the following picture for the composite 
spectrum during epoch 3: a high frequency (> 1.4 GHz) spectrum with a spectral index of —0.7 and a low 
frequency (< 1.3 GHz) spectrum with an index of —0.85. Comparing the composite spectrum at epoch 1 with 
this putative composite spectrum at epoch 3, we find that the low frequency spectrum has steepened while the 
high frequency spectrum has retained its spectral index while undergoing a secular decrease. 




Figure 2. (a) Radio spectrum of emission from the nova remnant of GK Per at different epochs (data from |2|, (s) (tJ). The 
lines connecting the points at frequencies above 1.4 GHz are for a spectral index of —0.7 whereas the line coimecting epoch 
3 data is for a = —0.85 and the line connecting the low frequency points of epoch 1 is for a = —0.4. (b) A schematic 
showing the timescales of acceleration and synchrotron losses which can explain the observed spectrum in epoch 3. Thus, 
the electrons at low energy should have larger losses whereas at the high energies, the acceleration process should balance 
(or be slightly larger than) the losses. 



We now attempt to analyse the effects of possible acceleration, synchrotron losses and adiabatic losses suffered 
by the electrons that could explain such an evolution of the composite spectrum. At high frequencies, if the 
acceleration timescale ta is close to (or slightly smaller than) the synchrotron loss timescale ts the spectrum 
will not undergo changes in its spectral index and will merely shift downwards (i.e., undergo a secular decrease) 
as a result of adiabatic expansion. At low frequencies, however, if ts < ta, it would lead to a steepening of 
the spectrum over time. Since the synchrotron loss timescale tg oc 1/E where E is the electron energy, the 
argument outlined above suggests a form for the acceleration timescale 

ta^E-", p<l (1) 

as shown in Fig|2jb). We denote the energy at which ta and tg intersect by Ec- The break frequency between 
the high and low energy spectra coiTesponds to the critical frequency of an electron with energy Ec- Fig|2jb) 
suggests that acceleration balances (or slightly dominates) over losses at electron energies > Ec while losses 
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dominate for electron energies < Ec- This could be a possible explanation for the evolution of the composite 
spectrum from epoch 1 to epoch 3. 

The time evolution of the electron distribution function / due to acceleration is commonly described by the 
following diffusion equation (e.g., Subramanian, Becker & Kazanas 1999) 



This equation describes the evolution of the distribution function as a result of diffusion in electron energy E. 
The acceleration timescale arising from such a mechanism can be written as (appendix A of |9J) 



If the energy dependence of the diffusion co-efficient is of the form V oc E^, the requirement on the accelera- 
tion timescale expressed in Eq ([l), when combined with Eq (|3} implies that 3 < ^. The physical mechanism 
accelerating the electrons must yield a diffusion co-efficient V which satisfies this requirement. 

4. Contribution to galactic cosmic rays 

The mininum energy in relativistic particles required to explain the observed synchrotron emission from GK 
Per using equipartition arguments comes out to be lO'^'^rj'^^'^ Joules where 77 is the ratio of energy contained in 
protons and electrons. This is fairly insignificant compared to the energy contained in particles in supernova 
remnants like Cas A (~ 2 x lO'^^rj'^^'^ Joules) which are significant contributors to the cosmic ray flux. Thus, 
although the nova rate of 30 yr~^ Js] is fairly high compared to the supernova rate of 0.02 yr~^ Q, the energy 
released per classical nova is insignificant and hence contributes little to the total cosmic ray flux. 
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